Magnetoelastic effects in multiferroic H0M11O3 
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We have investigated magnetoelastic effects in multiferroic HoMnOa below the antiferromagnetic 
phase transition TV w 75 K by neutron powder diffraction. The lattice parameter a of the hexag- 
onal unit cell of HoMnC>3 decrease in the usual way at lower temperatures and then shows abrupt 
contraction below Tjv whereas the lattice parameter c increases continuously with decreasing tem- 
perature and shows an anomalous increase below Tjv. The unit cell volume decreases continuously 
with decreasing temperature and undergoes abrupt contraction below Tat. By fitting the background 
thermal expansion for a nonmagnetic lattice with the Einstein- Griineisen equation we determined 
the lattice strains Aa, Ac and AV due to the magnetoelastic effects as a function of temperature. 
We have also determined the temperature variation of the ordered magnetic moment of Mn ion by 
fitting the measured Bragg intensities of the nuclear and magnetic reflections with the known crystal 
and magnetic structure models. 

PACS numbers: 61.05.fm, 65.40.De 



One of the most challenging subjects of the condensed 
matter physics is the coupling between spin and lattice 
degrees of freedom. The spin system in ferromagnetic 
or antiferromagnetic crystals is coupled to the ionic dis- 
placements via the dependence on distance of the ex- 
change interaction or spin-orbit or dipole-dipole inter- 
action. The coupling is called the magnetoelastic cou- 
pling and and resultant effects are called magnetoelastic 
effectPH^. The static part of this effect is the shift in the 
equilibrium ionic positions with resultant displacements 
of the phonon and magnon spectra. The dynamic part of 
this interaction is the magnon-phonon scattering and also 
their hybridization. The simplest static part of the in- 
teraction is the external magnetostriction, or the change 
in macroscopic crystal dimensions. This changes the lat- 
tice parameters below the magnetic phase transition and 
can be more easily determined experimentally. In ad- 
dition there are shifts in the atomic cooodinates within 
each unit cell. This is called internal magnetostriction 
which is more difficult to measure. The induced atomic 
displacements may or may not modify the symmetry of 
the lattice and also the anisotropic energy. There exist 
a large number of review articles that summarize the ex- 
perimental results and also to some extent some theoret- 
ical calculationPHS The magnetoelastic effect has drawn 
renewed interest in connection with the potentially use- 
ful colossal magnetoresistive and multiferroic materials. 
Here we have investigated the magnetostriction in mul- 
tiferroic hexagonal manganite HoMn03. 

HoMnOa belongs to the family of hexagonal mangan- 
ites RMn0 3 (R = Sc, Y, Er, Ho, Tm, Yb, Lu) that 
show multiferroic behaviouil^. These hexagonal man- 
ganites at paraelectric at high temperatures with the 
centrosymmetric space group PQ^/mmc. Below about 
1000 K they undergo a paraelectric-to-ferrielectric tran- 
sition to a non-centrosymmetric structure with the space 
group PQ^cm. At further lower temperatures of the or- 




FIG. 1: (Color online) Schematic representation of the hexag- 
onal crystal structure of RMnC>3 (R = Sc, Y, Er, Ho, Tm, 
Yb, Lu). Here we have used the structural parameters for 
HoMnC>3. The red trigonal bipyramids represent Mn ions 
surrounded by five O atoms and the blue spheres are Ho ions. 



der of about 100 K the magnetic hexagonal manganites 
order with a non-collinear antiferromagnetic structure 
with propagation vector k = 0. Among these hexagonal 
manganites YMnC>3 and HoMnC>3 have been investigated 
quite intensively^"^. We have performed neutron pow- 
der diffraction investigation of the temperature depen- 
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FIG. 2: (Color online) (a) Temperature variation of the 
diffraction intensities of HoMnC>3. (b) Temperature variation 
of a few low angle reflections from HoMnOs containing mag- 
netic contributions, (c) Temperature variation of the peak 
intensity of the 100, 101,102 and 104 reflections. 



FIG. 4: (Color online) Results of the refinement of the crystal 
and magnetic structures of HoMnOs at T = 2, 35 and 50 K. 
The upper tick marks below the diffraction pattern show the 
calculated positions of the nuclear peaks and the lower tick 
marks show the calculated positions of the magnetic peaks. 
The difference between the observed and calculated intensities 
has been plotted below the tick marks. 
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FIG. 3: (Color online) Results of the refinement of the hexag- 
onal crystal structure of HoMn0 3 at T = 100 K. The tick 
marks below the diffraction pattern show the calculated posi- 
tions of the diffraction peaks and the difference between the 
observed and calculated intensities has been plotted below the 
tick marks. 



dence of the crystal and magnetic structure of HoMnOs. 
Here we report the observation of considerable magne- 
toelastic effects below the Neel temperature T/v ~ 75 K. 

In order to determine quantitatively the magnetoelas- 
tic effect or the spontaneous magnetostriction one needs 
to measure very accurately the temperature variation of 
the lattice parameters and the unit cell volume in small 
temperature steps. These temperature variations often 
reveal anomalous behavior around the magnetic order- 
ing temperature. To extract the magnetoelastic effect 
one needs to know the background temperature varia- 
tion in the absence of magnetoelastic effect in a fictitious 
non-magnetic solid which otherwise resembles the com- 
pound under study. This is not an easy task and the 
disagreements between experimental results often arise 
from the procedure of determination of the background 
temperature variation of the lattice parameters. One 
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FIG. 5: (Color online) Temperature variation of the lattice 
parameters a, c, and the unit cell volume V of HoMn03 plot- 
ted on the left panel. The red curves in these figures represent 
the lattice parameter and the unit cell volume obtained by fit- 
ting the high temperature data by a Einstein-Griineisen equa- 
tion and extrapolated to the low temperature. This should 
give the background for the non-magnetic lattice. On the 
right panel the temperature dependence of the lattice strains 
Aa and obtained by subtracting the non-magnetic back- 
ground have been plotted. 



FIG. 6: (Color online) Temperature variation of the | ratio of 
HoMnOs plotted on the upper panel. The red curves in these 
figures represent the - ratio obtained by fitting the high tem- 
perature data by the Einstein-Griineisen function and extrap- 
olated to the low temperature. It is assumed that this gives 
the temperature variation of - for the non-magnetic lattice. 
On the bottom panel the temperature dependence of the - 
ratio obtained by subtracting the non-magnetic background 
have been plotted. 



of the commonly employed method is to fit the high 
temperature data above the ordering temperature by 
the Debye equation in the Griineisen approximation^ 
or the simpler Einstein equation in the same Griineisen 
approximation^ and then extrapolating to the lower 
temperatures to give the background variations in the 
absence of magnetism. The simpler Einstein-Griineisen 
equation^ serves as an excellent fit function for not only 
the volume but also for the individual lattice parameters 
and we have used it here. 

Neutron diffraction experiments were done on 
HoMnOs on the high intensity powder diffractometer 
D2C 19 of the Institute Laue-Langevin in Grenoble. The 
115 reflection from a Ge monochromator at a high take- 
off angle of 118° gave a neutron wavelength of 1.868 A. 



Approximately 5 g of HoMnOs powder sample was placed 
inside an 8 mm diameter vanadium can, which was fixed 
to the sample stick of a standard 4 He cryostat. We have 
measured the diffraction intensities from HoMnOs as a 
function of temperature in the range 2 — 300 K. Fig. [2] 
(a) shows the temperature variation of the diffraction in- 
tensities of HoMn03 . The temperature variation of a few 
low angle reflections from HoMnOs containing magnetic 
contributions is shown in Fig. [2] (b) and the temperature 
variation of the peak intensities of the 100, 101,102 and 
104 reflections is shown in Fig. 2] (c). The temperature 
variation of the intensities of the low angle reflections 
from HoMn03 clearly shows two magnetic phase transi- 
tions, one at the Neel temperature T/v ~ 75 K and the 
second spin reorientation transition at about T$r 40 
K. 
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FIG. 7: (Color online) (Upper panel) Temperature variation 
of the ordered magnetic moment m (Mnl) and m (Mn2) of 
Mn and m (Ho) of Ho ions of HoMnOs obtained by fitting the 
measured Bragg intensities of the nuclear and magnetic reflec- 
tions with the known crystal and magnetic structure models, 
m (Mnl) is the ordered magnetic moment corresponding to 
the spin rotated low temperature phase below Tsr ~ 40 K 
and m (Mn2) is the ordered magnetic moment of the higher 
temperature magnetic phase that develops below Tjv ~ 75 K 
and undergoes a phase transition below Tsr w 40 K to the 
spin reoriented low temperature phase. (Lower panel) Power- 
law fit of the temperature variation of m(Mn2) data close to 
T N . 



Rietveld refinemernpSI of crystal and magnetic struc- 
tures against the experimental diffraction data was done 
by the Fullprof program^. The refinement results from 
HoMnC>3 at T = 100 K in the paramagnetic state are 
shown in Fig. [3] At T = 100 K only the nuclear 
Bragg peaks were observed and the refinement of the 
nuclear structure was done by the know n crystal struc- 
ture of the hexagonal manganite d 12 l 13 l Fig. [4] shows 
the resul ts of refinement of the crystal and magnetic 
structureiPES at T = 2, 35 and 50 K in the three mag- 
netic phases of HoMnOa.The agreement factors and x 2 
values of these refinements are given in Table [TJ It is to 
be noted that magnetic structures of hexagonal mangan- 
ites including HoMnOa show homometry, which cannot 
be distinguished by powder neutron diffraction. Only po- 



larized neutron diffraction on single crystals can resolve 
the homometric structures^. 



TABLE I: Agreement factors for the refinement of crystal and 
magnetic structures of HoMnC>3 at different temperatures. R p 
and R wp are the unweighted and weighted agreement or R fac- 
tors whereas the R exp is the expected R factor corresponding 
to the statistics of the data and \ 2 nas the usual statistical 
meaning. 



T(K) R p (%) R wp (%) R exp {%) X A 
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8.13 


7.11 


3.15 


5.10 


35 


8.04 


7.03 


3.20 


4.82 


50 


8.53 


7.38 


3.22 


5.26 


100 


8.84 


7.35 


3.30 


4.98 



Figure [5] shows the temperature variation of the lat- 
tice parameters a, c, and the unit cell volume V of 
HoMnOs plotted on the left panel. The lattice parame- 
ter a decreases continuously with decreasing temperature 
in a very normal way, but close to the Neel temperature 
Tjv ~ 75 K shows the magnetoelastic or magnetostriction 
anomaly. The red curve represent the background vari- 
ation of the a lattice parameter for a nonmagnetic solid 
obtained by fitting the data in the paramagnetic state 
by the Einstein-Griineisen equation explained before. By 
subtracting the background from the data we determined 
the lattice strain Aa plotted on the right panel of Fig- 
ure [5] The c lattice parameter increases continuously 
on decreasing temperature down to Tjv and then shows 
abrupt anomalous increase or positive magnetostriction. 
By subtracting the background variation of the c lattice 
parameter (red continuous curve) from the data we de- 
termined the lattice strain Ac plotted on the right panel. 
Similar plots for the unit cell volume V and and the vol- 
ume strain AV are shown in Figure 5] in the left and 
the right panels, respectively. Fig. [6 shows in the up- 
per panel the temperature dependence of the - ratio. 
It varies linearly with temperature above about 100 K 
and below this temperature the ^ shows anomalous be- 
havior due the magnetoelastic effect. However we fitted 
the - ratio above 100 K by the Einstein-Griineisen equa- 
tion and extrapolated to the lower temperature. The 
fit is shown by the red line. Fig [6] shows in the lower 
panel the difference between the observed and the fit- 
ted ^ ratio. It is evident that the temperature variation 
of - also behave like the other lattice parameters with 
respect to the magnetoelastic effect. We note that al- 
though the magnetic ordering temperature of HoMnOs 
is Tjv ~ 75 K, the magnetoelastic effect starts becom- 
ing appreciable already below about T* m 100 K. This 
appearance of magnetoelastic effect already about 25 K 
above Tjv « 75K is due to short-range spin correlations. 
The spontaneous linear magnetostriction along the a-axis 
at T = is Aa/a — —2.45 x 10~ 4 and that along the 
c-axis is Ac/co = 1.31 x 10~ 4 . The spontaneous volume 
magnetorestriction in HoMnOa is AV/Vq = — 3.51 xl0~ 4 . 
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We can calculate the spontaneous volume magnetore- 
striction in HoMn0 3 from AV/V = 2Aa/a + Ac/c = 
-3.59 x 1(T 4 which very close to that obtained from the 
unit cell volume variation AV/Vq = —3.51 x 10~ 4 . This 
suggests that our method of extracting magnetoelastic 
effect is reliable or at least consistent. The magnitude of 
spontaneous magnetostriction in HoMnOa is relatively 
large and is comparable to those^l determined in the 
other multiferroic hexagonal manganite YMn03. 

Figure [7] (a) shows the temperature variation of the or- 
dered magnetic moments m of Mn ion corresponding to 
the high temperature magnetic phase of HoMn03 (shown 
by the red filled squares) and lower temperature spin re- 
oriented phase of HoMnOa (shown by filled blue circles) 
obtained by fitting the measured Bragg intensities of the 
nuclear and magnetic reflections with the known crys- 
tal and magnetic structure models. The ordered mag- 
netic moment of the high temperature magnetic phase 
decreases with increasing temperature above about the 
spin reorientation transition temperature T$r, ~ 40 K in 
the usual way and becomes zero at Tjy ~ 75 K. Below 
Tsr ~ 40 K the moment of this phase decreases and 
the magnetic moment of the spin reoriented low temper- 
ature phase increases abruptly and becomes saturated to 
about 3.1/zs at T = 2 K. The Ho magnetic moment is 
polarized at all temperatures below Tjv ~ 75 K but be- 
low about 50 K the Ho moment increases with decreas- 
ing temperature and becomes about 4/is at T = 2 K. 
The temperature variation of the ordered magnetic mo- 
ment of the Ho sublattice resembles the Brillouin function 
and indicates strongly the absence of any magnetic phase 
transition^ in H0M11O3 at lower temperatures. Figure 
[7] (b) shows the power-law fit of the magnetic moment 
of Mn in the higher temperature phase of HoMn03 be- 
low Tjy ~ 75 K. We have included the data points in the 
temperature range 50-74 K. The least squares fit of these 
data gave the power-law exponent j3 = 0.247 ± 0.005 and 



Tjv = 74.84 ± 0.05 K. Since there are no data close to 
T/v the the power-law exponent obtained cannot be as- 
sociated really with the critical exponent. However the 
determined Neel temperature Tjy = 74.84 ± 0.05 K from 
the fit is dependable within the experimental resolution. 
It is to be noted that the magnetoelastic effects shown 
in Fig. [5] in the lattice parameters are not sensitive to 
the spin-reorientation transition at Tsr ~ 40 K. It is not 
quite sure whether the ordered magnetic moment of the 
Ho sublattice at lower temperature influence the magne- 
toelastic effect substantially. 

The temperature variations of the lattice strains de- 
termined in HoMnOs do not look very smooth and do 
not become zero above Tjv ~ 75 K. This is partly due to 
the problem of determining the background temperature 
variation of the lattice parameters and unit cell volume 
by the Einstein-Griineisen function and also may be due 
to multiple magnetic phase transitions and the polariza- 
tion of the Ho magnetic moments. We could not therefore 
relate quantitatively the lattice strain and the magnetic 
order parameter as we did in YMnC<3 and other sim- 
ple antiferromagnetic systems^HUl We checked carefully 
whether we could get any information about the shifts in 
atomic coordinates in the unit cell below the magnetic or- 
dering temperature or the internal magnetoelastic effect. 
Although the agreement factors of the refinement of the 
crystal and magnetic structures against the data were 
reasonably good (see Table |l| , the resulting positional 
parameters do not show within the refinement accuracy 
any anomalous shifts in atomic coordinates below Tjv. 

In conclusion we have measured accurately neutron 
diffraction intensities of HoMnOs powder sample and de- 
termined the magnetostriction below the magnetic order- 
ing temperature. We have also determined the magnetic 
phase transitions and the temperature evolution of the 
magnetic moments of Mn and Ho ions. 
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